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I. INTRODUCTION 

Quenched disorder and dynamical disorder play an important role in the properties of many physical systems. 
Some examples that illustrate this role in a variety of very different contexts and support the wide interest in the 
subject of disorder include the propagation of fronts in porous media p], the conduction properties of doped semicon- 
ductors materials Q, the shift of the Curie temperature in ferromagnets 0, and the so-called Anderson localization 
transition 

At the same time, another broad field of perennial interest is pattern formation [f| because spatio-temporal struc- 
tures are omnipresent in the physical world |£g. Our own recent contributions in this field involve the discovery of 
novel mechanisms for the appearance of spatio-temporal structures upon periodic or random global alternation of 
pattern- free dynamics [7|. 

Herein we coadunate these two topics by investigating conditions under wich the presence of quenched disorder is 
the mechanism that triggers pattern formation. Thus our interest lies in a particular type of purely-noise-induced 
phenomenon: the appearance of patterns in disordered systems which in an ordered state exhibit no pattern formation. 
We concentrate on systems with quenched dichotomous disorder and coupling term a la Swift-Hohenberg Q , although 
one can easily envision generalizations of the formalism to other kinds of disorder and/or coupling terms. 

Noise-induced phenomena in spatially extended systems have been particularly active areas of investigation in the 
recent past 9] . Among these are phenomena involving pattern formation induced purely by fluctuations [ToL ITU IT2| . 
The word "purely" here emphasises the fact that the control parameter that determines the presence or absence of 
patterns is the noise intensity. Moreover, if the noise term is replaced by a non-fluctuating parameter, no spatio- 
temporal structures develop for any value of that parameter. Our study differs from these previous ones in that we 
consider quenched spatial disorder, that is, the fluctuations have no temporal dependence. 

Instead of focusing on a specific system at the outset, we explore some general conditions under which these systems 
exhibit purely disorder induced patterns. We then illustrate our findings with a family of systems that includes the 
paradigmatic models of noise-induced phase transitions 0,0,0] and noise-induced spatial patterns pTilTlif. We also 
identify a phenomenon in our system that has previously only been identified in systems with colored noise [l J. ITU IT^ | 
or field-dependent kinetic coefficients [T^. namely, a re-entrant behavior with increasing coupling strength. In other 
words, we find that increasing the coupling strength leads to non-monotonic behavior such that the patterns are most 
prominent for a finite value of the coupling and disappear altogether when coupling is too strong (or too weak). As in 
second and first order phase transitions in equilibrium systems, the re-entrance phenomenon can be either continuous 
(second order behavior) or it may present multi-stability and associated hysteresis (first order behavior). 

The paper is organized as follows. We introduce the formalism in Sec.[H] In Sec. IHII and Appendix lAl we present a 
generalized mean field approximation and state the requirements for pattern formation induced purely by disorder. The 
possible behaviors that can be deduced from these requirements are explored in Sees. HVl and Ivl Particular examples 
of systems that exhibit pattern formation are given in Sec. I VII and the order parameters used to characterize the 
patterns are introduced and related to one another in Appendix^! Numerical simulations that confirm the qualitative 
validity of the theoretical results are presented in Sec. IVIII Finally, we summarize the main results in Sec. IVIIII 

II. THE MODEL 

We consider the following stochastic dynamics for a scalar field </> r = </>(r, t) in the presence of dichotomous disorder: 
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Here £ r is a space dependent quenched dichotomous variable that models spatial disorder. The probability density of 
£ reads: 



p(0=P+*(£-A)+p_<J(£ + A), 



(2) 



that is, at any given site r the variable £ r takes on either the value +A or the value —A with probabilities p+ and 
P- respectively. The term C stands for the Swift-Hohcnberg coupling operator 



C = -D(kl + V 2 ) 2 - 
The effect of this coupling can be deduced by applying £ to a plane wave e ik r , 

Ce ikr = uj(k)e ikr , 



(3) 



(4) 



where u>(k) = —D(k^ — k 2 ) 2 is the continuous dispersion relation (we use bold for vectorial quantities and italic for 
their magnitudes). 

In order to implement a mean field theory for this system we need to distinguish from one another neighboring 
locations r and r', which in turn requires that we discretize the system. Since numerical simulations also involve 
discretization, this procedure does not interfere with the comparisons of theoretical and numerical results. With the 
understanding of the action of the translation operator 



exp \^ x ^j f( x ) = f( x + Sx )> 
it is straightforward to deduce a discrete version of the Swift-Hohenberg coupling operator £, 

-i 2 



C=-D 



AtJ ^ \ 2 dx. 



(5) 



(6) 



where d stands for the spatial dimension, Ax for the lattice spacing, and indicates a partial derivative with 
respect to component i of the position vector r — [x\, X2, ■ ■ ■ ,Xi, . . . , Xj). The continuum delta function 8{r — r') is 
replaced in the usual way by a ratio that contains the Kronecker delta and the lattice spacing, 6 r ,r' /(Ax) d . As in the 
continuous case, the discrete dispersion relation can be obtained by applying the operator © to a plane wave e 
to obtain 



ik-r 



Lo{k) = -D 



2 d 



-, 2 



(7) 



Here fcj denotes component i of the wave vector k = (k±, k%, . . . , ki, . . . , kd). 

Note that uj(k) is nonpositive for any value of k in both the continuous and discrete cases, and that in the discrete 
case it depends not only on the magnitude but also on the direction of k. Of particular importance in our subsequent 
analysis are those modes for which u>(k) = 0. In the continuum these are the modes with k — ko, which are all those 
that lie on a continuous hypersurfacc of radius kg around the origin. In the discretized system the magnitudes k* of 
the least stable modes are shifted from ko and depend on direction, as can be seen by solving Eq. J7J). The longest 
vectors such that cu(k*) = lie along the cartesian directions in k space, e.g. (k* , 0, 0, . . . , 0) and have magnitude 



,* 2 . (k Ax 

max k — — — arcsm 

Ax V 2 



The shortest lie along a reciprocal space diagonal, e.g. -^(k* , k* , k* , . . . ,k*), and have magnitude 



min k* 



Ax 



koAx\ 
2Vd J 



(8) 



(9) 



In our analysis in d = 2 we take ko = 1 and Ax = 1. The difference between maxfc* and min k* is smaller than 
3%, the two values being ir/3 — 1.0472 and 7r/3.0737 = 1.0221. Therefore it is only a mild approximation to neglect 
the directional dependence of the solutions of uj(k*) = and focus on the magnitude, w(k*) — 0. Furthermore, in 
simulations one must use a finite system of N d sites [i.e., of volume (NAx) d ], so that the allowed modes themselves 
form a discrete set, with each component separated from the next one by an interval 5k = 2ir/NAx. One way to pick 



3 



the least stable modes is to construct a ring of radius (k*) (which we shall simply call k* from here on) of thickness 
5k and to consider all the modes that lie in this ring. We can then estimate the number of modes n(k*) that are least 
stable by calculating the number of cells of volume (2ir/N) d in the ring: 

dn d / 2 fNAxk^Y' 1 , s 

"(fe*) = m I TT I nZ ) ■ (10) 



r(d/2 + l) V 2tt 

Although slight variations in the particular way of counting are possible, for sufficiently large N the differences are 
small 



III. MODULATED MEAN-FIELD THEORY 



To establish the existence of patterns of a characteristic length scale in the steady state, we seek a spatially periodic 
structure defined by wavevectors whose magnitude is associated with the inverse of this length scale. The appropriate 
wavevectors to focus on are precisely those of magnitude k* , that is, those for which ui(k) = 0. A mean field theory 
requires that we make an ansatz about the behavior of the stationary field at sites r' ^ r which are coupled to the 
focus site r by the operator C We require that the ansatz capture correct limiting behaviors and also incorporate an 
appropriate spatial modulation. Our choice is 



A(k*) 2^ cos [k- (r-r')] 
{fe«} 



B, 



(11) 



where the sum (or, in an infinite system, the integral) is over wavevectors of magnitude k* . The constant B is 
specified below. Our ansatz incorporates the assumption that all modes of magnitude k* contribute with equal 
(direction-independent) weight A(k*). In Appendix lAl we show that the action of the coupling operator on this ansatz 
is given by 



£<p r = Di [n(k*)A(k*) - r ] + B(D 1 - Dhfi). 



where 



D\= D 



2d 
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(12) 



(13) 



Substitution of Eq. i|12|) in Eq. QJ then leads to an equation that depends only on a generic site index r that can 
simply be dropped: 



4> = f{4>) + 9m + Di [n(k*)A(k*) -</>] + B(D X - Dk 4 ). 
In the steady state we can write the explicit equations associated with Eq. (|14fl as 

= F + (4>) + Di[n(k*)A(k*)-(f>] + B(D 1 -Dk^), 



(14) 



= F_(4>)+D 1 [n(k*)A(k*)-4>]+B(D 1 -Dk%), (15) 

where we have introduced the shorthand notation 

F ± (4>)=f^)±g((j>)A, (16) 

We denote the solutions of Eqs. I|15|) by <f>± respectively. 

The amplitude A(k*) and the constant B are the mean field quantities that must be chosen self-consistently to 
complete the solution of the problem. To close these equations we choose 

B= P+ 4> + +p^ (17) 

where the <f>± are the steady state fields for each of the separate dynamics in completely ordered systems, that is, one 
with p + = 1 and one with p_ = 1 (see below) . Furthermore, we impose the self-consistency condition 

/oo 
d<t)4> P (4>;A(k*))d4>. (18) 
-oo 
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Since 

p{fr A{k*)) = p+5(<t> - 0+) +p_<5(0 - </,_), (19) 

the self-consistency condition can be rewritten as 

n(k*)A(k*) - $+) +p-(<t>- - ?_), (20) 

where it should be stressed that the <j>± are of course functions of A(k*). Finally, Eqs. I|15(l can then be rewritten as 

= F±fo±) + D iPt (0 t - <t>±) - Dkt{ P+ 4>+ + P -4>-). (21) 

We are particularly interested in systems in which there are no patterns in any ordered state, that is, where patterns 
are purely disorder-induced (we will choose the functions / and g accordingly) . For each value of the dichotomous 
disorder parameter we can write a deterministic evolution equation for the dynamics. We insist that each of these 
evolution equations describe a system that intrinsically has at least one steady state and hence we insist that the 
associated forces be confining. We also require that each of these dynamics be pattern-free, so that there must be 
exactly one steady state solution for each, and this steady state solution in each case must be a constant independent 
of r. The confining condition requires that for all A 

lim F± (<p) - Dk£<f» = oo 

lim F ± {4>) - Dk^<f) = -oo. (22) 

Indeed, it is reasonable to require that the confinement not be due simply to the coupling, and thus to require that 

lim F±(<f)) = oo, lim F±(<j>) = -oo. (23) 

<p — v—oo 4> — >-oo 

The steady state conditions follow directly from Eq. (JTJ: 

= F ± (<P) - DkU- (24) 

The solutions [already introduced in the choice of B in Eq. Q17[l] are denoted by 0+ and 0_ respectively. Further, to 
avoid any insertion of patterns other than those induced by disorder, we also insist that f(<p) itself be associated with 
only a single steady state, i.e., that the equations /(</>) — Dk^cj) = and /(</>) = also have only a single solution. On 
the other hand, if we approximate the dynamics in these ordered systems by our mean field ansatz, and insist that 
A(k*) — since there are no patterns, the equations in the steady state with B chosen as in Eq. (|17fl in each case 
read 

U = F±{4>)-D l ct> + 4> ± {D 1 -Dkt). (25) 

Clearly, the solutions of the two equations are again <j) + and <f>- respectively, and thus our choice of B is consistent 
with the exact unpatterned solutions in the ordered steady states. 

For simplicity, in most of the remainder of this work we set p+ = p- = 1/2, although we hasten to add that this 
condition is not necessary for the appearance of patterns. To support this statement, we show in Fig. nt ne order 
parameter S(k*) = n(k*)A 2 (k*) (discussed in more detail in Sec. I VII and Appendix lB|l vs p+ for a particular model 
considered later; the details of that model are not important at this point. A nonzero value of the order parameter 
indicates the appearance of patterns, and a higher value is indicative of more pronounced patterns. In this particular 
instance the strongest patterns occur whenp + = 1/2, but the figure shows that other values of p + also lead to pattern 
formation for the same model and parameter values. Furthermore, we also impose the requirement that /(</>) be an 
odd function and g(4>) an even function of <f). This immediately leads to the symmetry 

F + (0 + ) = -F_{-4>_) (26) 

and, associated with the pattern-free solution A(k*) = 0, 

4>+ = -0- (27) 

In particular, the average homogeneous solution B = 0. This is always one of the possible solutions of the problem. 
The stability of this state, and the possible appearance of other solutions with A(k*) ^ 0, are the subjects of our 
further analysis. 
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FIG. 1: Order parameter vs p+ for a model detailed in Sec. I VII Pattern formation is associated with a nonzero value of the 
order parameter. 

We collect, then, the system of equations to be analyzed with these simplifications. Furthermore, taking advantage 
of the symmetries established so far, the equations take on a more tractable appearance if we define new variables 

x = —4>+, y= — 2"0-, xq = —4>+^ — — (j>-. (28) 



Thus the equations to solve are [cf. Eq. Ijl5(l ] 



with 



F(x) — x = y 

F{y)~y = x (29) 
B^{k*)A{k*) =x-y, (30) 



where we have applied the symmetry F+(x) = —F_(y) and have dropped the subscript + since we only need to use 
F+: 

F(x) = F + (x)=f(x)+g(x)A. (31) 
To this we add the requirement that follows from Eq. (|2*3^l . 

lim F'(x) < 0, (32) 

x — *±oo 

where the prime denotes a derivative with respect to the argument. The solution x = y = xq is pattern- free. We seek 
solutions with x ^ y to establish the appearance of patterns. Note that if a pair (x,y) solves Eqs. <|29[) . then so does 
the pair (y, x), simply leading to a reversal in the sign of A(k*). Since negative values of A(k*) can be interpreted in 
terms of an overall spatial phase, the second pair adds no new physical information beyond the symmetry statement. 

The stability of the pattern-free solution becomes an important issue in our further discussion. This can be 
established by a linear stability analysis around x = y = xq. The time evolution of small perturbations about this 
solution is obtained by expanding the evolution equations 

F(x) -x-y = x, F(y) - y-x = y, (33) 

at (x, y) = (xq + Sx, xq + Sy) and retaining terms up to first order in the perturbations: 

\SyJ \ -1 F'(x ) -l)\8y) 
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The 2x2 evolution matrix has eigenvalues 

A+ = F'(x ), X- = F'(x Q ) - 2. (35) 

The largest eingenvalue (largest Lyapunov exponent) is clearly A+, and it determines the stability of the pattern- free 
solution: 

F'(xq) < solution is stable 

(36) 

-F'(xo) > solution is unstable. 

IV. PATTERN FORMATION INDUCED PURELY BY DISORDER 

We seek solutions for Eqs. I|29(l and, in particular, solutions with x ^ y. For given potential functions, one could 
produce three-dimensional plots of F{x) — x — y and F(y) — y — x vs x and y and observe the intersections of these 
surfaces. A more intuitive graphical way to organize this search is presented in the three panels of Fig. EI The origin 
represents the pattern-free solution, and this is the only point on the line y = x that solves the equations since F(x) 
is not an odd function. Since Eqs. I|29ll are invariant under the transformation x < — > y, the line x = y defines a 
specular plane, that is, F(y) — y is a specular image of F(x) — x with respect to that symmetry plane. The asymptotic 
behavior I1M2I) tells us that for sufficiently large \x\ in the upper left quadrant, F{x) — x lies above the line x = — y 
and has a slope < — 1. This is schematically indicated by the thick solid line in the upper left quadrant of each panel 
in the figure. Similarly, in the lower right quadrant the thick solid line recognizes that F(x) — x must lie below the 
line x = — y with a slope < — 1. The specular symmetry around the line x = y then leads us to the asymptotic thick 
dashed lines representing the behavior of F(y) — y. 

In panel (a) we illustrate a case for which the pattern- free solution is stable, that is, F'(xo) < 0. This is indicated 
by the thick solid line going through the origin. Again, the specular symmetry leads us to the thick dashed line to 
indicate the approriate slope for F(y) — y. Now it is clear that the thick solid lines can be connected, and the thick 
dashed lines can be connected, in such a way that the two lines do not cross anywhere else but at the origin. Thus, 
when the state (xq,Xq) is stable, there may not appear any other stationary states and the system may simply be 
pattern free. This is the case we have sketched in panel (a). 

On the other hand, it is possible to connect the solid and dashed curves respectively in such a way that there 
are two additional crossings of the curves (actually four crossings, but only two provide independent information). 
These represent two additional steady state solutions, each leading to pattern formation since x ^ y. However, only 
one of the two is stable. This is illustrated in panel (b). This represents the case of coexisting stable states, one 
pattern-free and the other patterned, separated by an unstable state. Such coexisting states are characteristic of 
first order phase transitions. One can carry this further and envision further crossings, always an even number of 
additional crossings, representing stable patterned states that coexist with one another and with the pattern-free 
state, separated by unstable states. 

In panel (c) we illustrate a case for which the pattern-free solution is unstable, that is, F'(xo) < 0. Again, this 
is indicated by the thick solid line going through the origin together with its dashed specular partner. Now it is 
clear that a connection of the lines necessarily leads to at least one crossing with x ^ y. In other words, when the 
pattern-free state becomes unstable, at least one patterned state is necessarily stable. Therefore, a sufficient condition 
for the occurrence of patterns is that the pattern-free state become unstable. 

In summary, we have found the following general behavior: 

• When the pattern-free solution is stable, there may or may not occur one (or more) stable solution(s) that leads 
to pattern formation, the patterned and unpatterned stable states being separated from one another by unstable 
solutions. The appearance of such an additional stable solution is indicative of a first order phase transition, 
with the usual coexistence and hysteresis characteristics. 

• On the other hand, when the pattern-free solution is unstable, a patterned stable state necessarily appears. The 
destabilization of the pattern-free state might simply mark the end of the coexistence region of a first-order 
transition described above, or it might mark the occurrence of a second order transition. These alternatives are 
discussed in further detail below. 

We note here that although we have not said so explicitly, it is implicit in this entire analysis that patterned states 
can only exist if the potential functions are nonlinear. Linear forms can not satisfy the conditions that lead to the 
emergence of patterns. 



7 




FIG. 2: Schematic of possible solutions of the mean field equations. The origin represents the pattern-free solution. Other 
intersections of the thick solid curve (which represents F(x) — x) and the thick dashed curve (which represents F(y) — y) are 
solutions of Eq. 12' H that lead to patterns. Panel (a) represents a case in which the pattern-free solution is stable and there are 
no other solutions, that is, there is no pattern formation. Panel (b) describes the coexistence of the stable pattern-free state 
with a patterned stable state. In panel (c) the pattern-free state is unstable and the only stable state is a patterned state. 



V. MEAN FIELD SOLUTION - GENERAL FEATURES 



Another way to exhibit the variety of possible transitions discussed above is to focus on the explicit solution of 
the mean field equations and consider the resulting values of A(k*). We start with the first Eq. (|2Tj)) and subtract 
x from both sides, F{x) — 2x = y — x. Similarly, we subtract y from both sides of the second equation, to write 
F(y) — 2y = x — y. Using Eq. (|30|l then implies that we can write 

D in {k*)A(k*) = -F{x) + 2x = H-{x) 

= F(y)-2y = H + (y). (37) 

We can invert these relations, 

x = HZ 1 (A(k*)), y = H- 1 (A(k*)). (38) 

The self consistency condition (|30|l for the mean field solution can then be written as 

Dm(k*)A(k*) = [HZ^Aik*)) - Hl\A(k*))] = D!ii(k*)G(A(k*)). (39) 

Since A(k*) = is always a solution (albeit not always stable), we know that G(0) = 0. Furthermore, the symmetry 
of the problem implies that if A(k*) is a solution, then so is —A(k*). Therefore G(A(k*)) is an odd function. An 
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FIG. 3: Schematic of G(A) vs .4. Intersections of the thick solid curve and the thick dashed line represent solutions of the 
mean field equations. Each panel represents the same situation as the corresponding panel in Fig. Panel (a) depicts the case 
of a stable pattern-free solution and no other. Panel (b) shows the case of coexistence of the stable pattern-free state and a 
patterned stable state (separated by an unstable solution). In panel (c) the pattern- free state is unstable and only a patterned 
state is stable. 



expansion about zero thus has only odd powers, 

G(A(k*)) = aA(k*) + bA 3 (k*) + O (A 5 (k*)) 



(40) 



Our analysis proceeds on the basis of the first term of this expansion as well as the asymptotic behavior of the 
function G. The graphical representation of this analysis is shown in the three panels of Fig. |3 each associated with 
the corresponding panels in Fig. [5] 

Let us first deduce the asymptotic behavior of G by considering the slope 



dHZ\A) 
dA 



1 



1 



d.r 



H_{x) 



jL 

dx 



[-F(x) + 2x] I 2 - F'fc) 



where we have applied the general relation 



—h~ l ( \ = 1 

dz ' {Z > h'[h- l {z)] 



(41) 



(42) 
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and the prime, as usual, denotes a derivative with respect to the argument. In particular, we thus find that 

lim — - — - = lim — — - < — , (43) 

^±00 dA a:-±oo[2-F'(s)] 2' V ' 

where we have used Eq. (|32(l . The slope — dH+ 1 (A)/dA leads to exactly the same asymptotic result, so that it follows 
from Eq. that 

.. dG(A(k*)) 1 fAAS 
lim a \. < 1. (44) 

.A—±oo dA{k*) y ' 

This clearly implies that asymptotically the function G(A) in the positive half-plane must lie below the diagonal line, 
as we have drawn in the three panels in Fig. [31 (in the negative half plane it must lie above, again as shown). We keep 
in mind that patterned solutions are associated with intersections of the function G(A) and the diagonal, away from 
the origin. 

Next we look at the behavior of G near the origin, and consider the coefficient a in the Taylor expansion 140|l . 
The calculation of this coefficient involves precisely the steps followed above, but with the functions evaluated at xq 
instead of asymptotically. We readily obtain 



dG(A(k*) 



dA(k* 



Am=o 2-F'( Xo y (45) 



Now suppose that F'(xo) < 0, a condition that according to Eq. Ij36(l means that the pattern-free solution is stable. 
The slope of G(A(k*)) near the origin is below the diagonal, a < 1. It is then possible that other than the crossing at 
the origin there is no other crossing, that is, only the stable pattern-free solution exists. This is illustrated in panel 
(a) of Fig. |3 and corresponds to the situation in panel (a) of Fig. |2 

Another possibility is that there are an even number of additional crossings, as illustrated in panel (b). One of the 
two additional solutions shown in the panel would be stable, the other unstable, and the stable patterned solution 
would coexist with the pattern-free solution under the circumstances shown in the panel. Again, this corresponds 
to the situation in panel (b) of Fig. [21 Below we establish further conditions on the forces that might lead to this 
behavior. 

Next suppose that F'(xq) > 0, the condition that according to Eq. (|36|l is associated with an unstable pattern- free 
solution. The slope of G(A(k*)) near the origin is now above the diagonal in the positive half-plane, a > 1, and 
another crossing than the one at the origin certainly occurs, thus insuring a stable patterned solution. This is shown 
in panel (c) of Fig. [31 and corresponds to panel (c) in Fig. El 

Although neither a necessary nor a sufficient condition, it is apparent that the additional crossings in panel (b) of 
Fig.[31might be accompanied by a positive curvature of G{A) near the origin, as shown in the panel, that is, additional 
crossings might occur if the coefficient b in Eq. I|4(J|I is positive. Likewise, case (c) is likely to be associated with a 
negative curvature near the origin, b < 0. 

It is useful to take cognizance of the possible sequences of behavior as one varies model parameters. For instance, a 
sequence (a) — > (b) — > (c) would signal a first order phase transition from an unpatterned state through a coexistance 
regime of unpatterned and patterned states to a regime where only the patterned state is stable. On the other hand, 
a sequence (a) — > (c) would represent a second order phase transition from an unpatterned to a patterned state. 
Transitions might be re-entrant, so that a returns (c) — > (b) — > (a) (first order) or (c) — * (a) (second order) are 
possible. In the next section we explore these results in the context of particular examples. 

VI. MEAN FIELD SOLUTION - PARTICULAR EXAMPLES AND PHASE DIAGRAMS 

Consider the particular family of force functions 

/(<« = -0(1 + ^ 2 r, 5 w.) = (i + <m (46) 

with m > 1. We will concentrate in particular on the cases m — 1 and m — 2. 

It is straightfoward to check that with this family of force functions, for any value of A and coupling D, Eqs. i|24|) 
have only a single real solution, that is, each of the ordered systems is pattern-free in the steady state. Thus, any 
pattern observed in the disordered system is purely a consequence of the disorder. 

We start with m — 1, a choice that has been made in a number of studies of purely noise induced transitions p^.ll4l 
Il5| . This is a particularly useful example because it can be solved analytically in mean field. Furthermore, We solve 
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FIG. 4: Phase diagram for the case m = 1. The shaded part of the diagram is a coexistence region and the unshaded interior 
denotes the occurrence of a single patterned phase. See text for a detailed description. The inset shows the behavior of the 
order parameter defined in the text as a function of D along the line A = 2. 



the coupled equations (|29|l fully and calculate the resulting amplitude A using Eq. (|30[1 , check the solutions for their 
stability properties, and thus construct a phase diagram as a function of the parameters A and D. This procedure 
yields the diagram shown in Fig. 0] We will describe the features of the diagram in terms of the analytic solutions to 
the problem as well as the general diagnostic measures described in the previous sections. 

Equations l|29l) are cubic and yield altogether five solutions. One is the pattern-free solution xq = yo f° r which 
A(k*) — 0. Of the remaining four solutions, only two are distinct (the other two are their negatives) and they yield: 

n(k*)A ± (k*) = ±- (A 2 ± A^2(D 1 - 2) - ±D t - lj . (47) 

Note that in our equations we continue to use D\ for aesthetic reasons, whereas our illustrations involve D. In two 
dimensions with the parameters specified earlier, the two are related by D\ = 13D. Several features of these solutions 
are noteworthy and describe the results in Fig. 

1. Both solutions are complex if D\ < 2. Therefore, the only stationary state when D\ < 2 is pattern-free (panels 
(a) in Figs. |3 and EJ). 

2. When D\ > 2, both solutions are real in the parameter range 

A > i (V2OD1 - 2) + . (48) 

One of these solutions (-4+) is stable, the other unstable, and this regime marks the shaded wedge in the figure, where 
the stable patterned solution coexists with the pattern-free solution. Entry into this region marks a first order phase 
transition (panels (b) in Figs. [21 and OJ). The two solutions merge (and the wedge closes) at the point A = 2, D\ = 2. 

3. When D 1 > 2 and 

i (-V2(£>i - 2) + 2^2^) < A < i (V2(£>i - 2) + 2^2^) , (49) 

only the solution A+ is real. This solution is stable, and delimits the unshaded region within which there is only a 
single steady state, which is patterned. The transition into this regime is part of the first order phase transition if the 
crossing is from the shaded region, or a second order transition if from the pattern-free region (panels (c) in Figs. |2 
and|3J). The boundary of this regime is thus precisely the curve defined by the bounds given in Eq. I|49() . Alternatively, 
we can invert this equation and express these same bounds as 



£>i = I (5A 2 - 3 ± 4A VA 2 - 3) 



(50) 
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FIG. 5: Phase diagram for the case m — 2. As in the case m — 1, the shaded portions represent coexistence regions and the 
unshaded interior indicates a single patterned phase. The inset graph shows the order parameter as a function of D for A = 2. 
Note the double-hysteresis behavior as a function of the coupling. 



Note that the two curves meet at the point A = -\/3, D\ = 8/3, which is the leftmost point of the contour. 

4. According to this description, the boundary of the unshaded region is precisely the curve along which F'(xq) = 0, 
thus bounding the regime within which the pattern-free stationary state becomes unstable. This is indeed the case, 
as we can establish without solving explicitly for xo as a function of D\ and A (the result is rather cumbersome). 
In terms of our original variable, for our potential we have the explicit expression F'(<fr) = —3(j) 2 + 2A(f> + 1. Setting 
this to zero yields 4>+ — (A ± ^A 2 — l)/3. Requiring that these values indeed be ones that solve Eqs. I|2HJ1 with 
x = y = xq or, in our original notation, F(<p + ) — D\(j) + , again gives exactly Eq. I|50[) . 

We call special attention to the striking re-entrance behavior of pattern formation as a function of the coupling 
D: sufficiently strong coupling destroys any patterns. Note that this implies that for a given value of the disorder 
parameter A there is an optimal coupling for which the patterned structure is most pronounced. A more nuanced 
discussion of this behavior requires quantification in terms of order parameters. In Appendix [B] we introduce in 
Eqs. HB6|) and l|B7|l the total power spectrum S(k) and the flux of convective heat J. They are related to one another 
in Eq. I|B8I) . and they both contain useful information. In general, when there is no pattern at all, >S'(A;) is independent 
of k and of 0(J/N), i.e., S(k)/J = 0(1/N) for all k. On the other hand, when there is a pattern of wavevector 
magnitude k* , then S(k*) is of O(J), i.e., S(k*)/J — 0(1). A larger value of S(k*) indicates a stronger pattern. On 
the other hand, a larger ratio S(k*)/J indicates a more coherent pattern. It is possible (as we shall see below) for a 
pattern to become more (or less) coherent even as it becomes less (or more) pronounced. 

In our mean field theory, however, we do not have access to all of this information because we do not deal with all 
modes but only with those that dominate the pattern when one is present. We are therefore restricted to choose as 
an order parameter the value of the power spectrum at k*: 

S(k*)=n(k*)A 2 {k*). (51) 

Within our theory S (k*) and J are the same. When there is no pattern they are both zero and their ratio is undefined. 
When there is a pattern, S(k*)/J = 1. Our theory can therefore predict whether a pattern will become stronger as 
parameters are modified, but not whether it will become more or less coherent. 

The inset in Fig. 0] shows the power spectrum for the m = 1 model as a function of coupling at A = 2 obtained 
from our theory. The values are just at the edge of the second order phase transition region and, accordingly, the 
order parameter rises continuously (albeit steeply) from zero upon entry into the patterned region. It also vanishes 
continuously as we exit the patterned region, again indicating a second order transition from the patterned state back 
to a pattern-free state (re-entrance) with increasing coupling. The most pronounced pattern is associated with the 
maximum in this curve, which can be found from Eq. (|47|l to occur at D\ = 2.94. 

We have also carried out the same set of calculations for the potential functions H46fl with m = 2. While some 
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phase boundary information can be obtained analytically, full analytic solution is no longer possible because the 
equations are now quintic. A numerical solution is straightforward and leads to the phase diagram shown in Fig. [S] 
The inset shows the order parameter for A = 2, where the appearance and disappearance of patterns with increasing 
coupling are both first-order transitions. As a result, the order parameter curve jumps discontinuously at both ends 
and exhibits, at both ends, the usual hysteresis behavior associated with first-order transitions. 

VII. NUMERICAL SIMULATIONS 

In order to check the predictions of the modified mean field theory we perform numerical simulations for the family 
of force functions given in Eq. (|46[) with m — 1 and m = 2. We implement an Euler numerical scheme on a two 
dimensional square lattice with periodic boundary conditions. The values of the parameters used in the simulations 
are Ax = 1, fco = 1, L = NAx — 64, and At = 10~ 3 . The aspect ratio, which measures the number of wavelengths 
of the expected patterns that fit into the system, is L/ (27r/fc*) ~ 10. Below we typically present averages over ten 
disorder configurations. 

Pattern formation is characterized by the total power spectrum, S(k*) [cf. Eq. jB6j)], where the sum runs over 
wavevectors k whose magnitude k lies in the interval [k* — 2ir/L,k* + 2Tr/L}. In addition to this parameter and the flux 
of convective heat J, another useful quantity for characterizing the system is the angular average of x s ) = <pk<P-k- 

The first panel in Fig. [S] shows the order parameter S(k*) as a function of the disorder intensity parameter A for 
the family m — 1 and for coupling coefficient D = 5. In the absence of disorder there is no pattern, S(k*) as 0, but a 
pattern clearly develops as A increases and is therefore entirely due to the disorder. The inset graphs show x (k) by 
means of a wavevector density plot for A = 3.5 and for A = 5. Clearly, a ring of unstable modes develops around k* , 
and the ring becomes more prominent with increasing disorder (as measured by the value of A). While we recognize 
that our modified mean field theory does not predict the transition values quantitatively (for D = 5 we predict a 
patterned state to first appear when A = 11.5 while the simulations already show pattern formation for much smaller 
values of A), the qualitative behavior is as predicted. It should be noted that while the intensity of x(k) at the most 
unstable modes increases as indicated by the gray scale, the width of the ring around k* also increases with increasing 
A. While a higher intensity indicates a more pronounced pattern, the width is a measure of the coherence of the 
patterns, increasing width indicating greater decoherence. As mentioned earlier, the ratio S(k*)/J can be used to 
characterize the coherence of the patterns. This nonmonotonic ratio is shown in the second panel in Fig. As also 
noted earlier, this ratio can not be obtained from our mean field theory. 

The re-entrant behavior as a function of the coupling predicted by the modified mean field theory is shown in the 
first panel of Fig. [7| for m — 1 and A = 2.5. The inset panels depict density plots of x(k) for D = 0, D — 2, and 
D = 12. Again, the agreement with the theory is only qualitative but nevertheless dramatic because re-entrance with 
increasing coupling strength is a rare phenomenon. Notice the destabilization of a ring of modes around k* with 
increasing coupling, and its subsequent intensity fade-out as D increases further. It should be noted that not only 
does the intensity of x(k) at the most unstable modes decrease as indicated by the gray scale, but the width of the 
ring also shrinks around k* as the coupling grows. The patterns thus become less pronounced but more coherent with 
increasing coupling. A representation of S (k*) /J, shown in the second panel of Fig. [71 reveals that although the 
system presents a re-entrant behavior with the coupling, the coherence actually increases monotonically as a function 
of D. 

The actual stationary spatial patterned configuration induced by the disorder is shown in the lower panel of Fig. |S| 
for m = 1, A = 2.5 and D = 5, i.e., the rightmost point in the first panel of Fig.[7| The upper panel in Fig. shows 
the particular configuration of quenched disorder in the simulation that leads to the stationary field shown in the 
lower panel. 

We have thus confirmed the disordered-induced pattern formation phenomenon but have not yet ascertained our 
other prediction, namely, the occurrence of hysteresis in some cases. We predict hysteresis to occur in the shaded 
regions of the phase diagrams in Figs. 01 and [SJ Although we do not necessarily expect quantitative agreement with 
the particular values of parameters that lead to hysteresis, the question is whether hysteresis is observed at all. 

Hysteresis implies a memory of the initial conditions. To detect hysteresis we carry out two simulations. In one, 
starting from an unpatterned initial condition (f> r (to) = for all r, we calculate the total power spectrum in the 
steady state as a function of increasing coupling. Then we decrease the coupling along the same phase space path, 
but now our initial condition for each value of the coupling is the stationary state obtained in the simulation with 
the previous value of D. A difference in the spectrum obtained for these two different conditions reflects hysteresis 
and the attendant coexistence of different stationary states (one patterned and one unpatterned, as discussed earlier). 
In Fig. we observe precisely this behavior for the family of functions m = 2 when increasing (black circles) and 
decreasing (white circles) the coupling for a fixed value of the intensity of the disorder, A = 3. Hysteresis is observed 
between D = and D ~ 7. The qualitative behavior is thus exactly as predicted by the mean field model. 
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FIG. 6: First panel: S(k*) vs A for m = 1 and D = 5. The insets are density plots of x(k) for A = 3.5 and for A = 5. Note 
the destabilization of the fe*-modes as A grows. The gray scale used in the density plot is the same for both insets. Second 
panel: S(k*)/J vs A. Note the nonmonotonic behavior, indicating maximal coherence at around A = 2. 

The coexistence of patterned and unpatterned states is discernible in the density plot insets of \ s ) m the figure in 
that they are no longer cleanly annular but now include contributions from wavevectors other than those of magnitude 
near k* . The difference between the two insets (indicating different relative contributions of patterned and unpatterned 
states) is visible but, an even clearer rendition of the difference is seen in the angular average (x(k)) e shown in the 
other two inset panels. Note, for instance, that the contribution of modes near k — are relevant when increasing 
the coupling but very small when decreasing D. 



VIII. SUMMARY 



Using a modified mean field theory, we have explored general conditions under which spatially extended systems with 
coupling a la Swift-Hohcnbcrg exhibit pattern formation purely induced by the presence of quenched dichotomous 
disorder. We have illustrated the phenomenology with a family of force functions that includes the paradigmatic 
models of noise-induced phase transitions and noise-induced spatial patterns, among them one that can fully be 
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FIG. 7: First panel: re-entrant behavior of the total power spectrum as a function of the coupling for m = 1 and A = 2.5. 
The inset panels show x(k) for D = 0, D = 2, and D — 12 by means of density plots: the same gray scale is used in all 
cases. Although S(k*) shows non-monotonic behavior, the pattern coherence increases as the coupling grows, as seen in the 
monotonic behavior of the ratio S(k*)/J in the second panel. 

solved analytically within the mean field model. We show that pattern formation can be achieved through continuous 
(second order) and discontinuous (first order) transitions, and that the pattern formation phenomenon is re-entrant 
with the coupling. Thus, increasing the coupling eventually destroys the order. All these predictions have been 
checked by means of numerical simulations, and we find full qualitative (albeit not quantitative) agreement with the 
theory. 

Beyond the capabilities of the modified mean field approach, we have explored not only the occurrence of patterns 
but also their coherence. The numerical simulations show that even while patterns become stronger (as reflected by a 
larger value of the total power spectrum at the particular wavevector magnitude characteristic of the pattern), they 
may become less coherent (as reflected by a larger contribution of neighboring wavevectors). The converse is also 
possible: patterns that become more coherent even as they weaken. 

We end by stressing that the methodology and formalism developed herein can be generalized straightforwardly to 
other kinds of disorder and/or couplings. 
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FIG. 8: Density plot showing a stationary pattern purely induced by disorder for m = 1, A = 2.5, and D = 12 (lower panel). 
The upper panel shows the underlying configuration of quenched disorder for this particular simulation. 
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FIG. 9: Total power spectrum as a function of the coupling for m = 2 and A = 3. The hysteresis is reflected in the difference 
between the black and white circle phase space trajectories, obtained when increasing and decreasing D respectively. The insets 
show density plots of x (k) and their angular averages. 
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APPENDIX A: GENERALIZATION OF MODIFIED MEAN FIELD THEORY 

We begin by illustrating some detailed dependences associated with the ansatz field l|ll|) and the action of C on it. 
For example, for an r' that is m lattice sites away from r = [xi, x 2 , ■ ■ ■ , Xd) in the direction j the ansatz reads 



(xi, x 2 , ■ ■ ■ , Xj + mAx, . . . , Xd) = A(k*) cos(mAxkj) + B. 

{k*} 



For an r' that is in the immediate positive diagonal location away from r we have 
<j)(x\ + Ax, X2 + Ax, . . . , Xj + Ax 



Next, to apply the discrete version JfJJ) of L we must elucidate the effect of the operators Y^t=i sm h 2 ' A " 



2 dx 



the field <p r for n = 1, 2. With n = 1, we use the relation 2 sinh (y/2) = [cosh(?/) — 1] and note that 

d 

i 1 1 \ i 

i X2 , • ■ • , Xj , . . . , Xd) 



f d \ 1 
^2 cosh ( Ax — J <j) r =- [0(xi + Ax, 



+ 4>{xi - Ax,x 2 , ■ ■ ■ ,Xj, . . . ,x d ) 



{xi,x 2 , ...,%• + Ax, . . . ,x d ) 



(xi,x 2 , ...,Xj- Ax, ...,x d ) 



(Al) 



. . . . , x d + Ax) = M k *) cos[Ax(fci + k 2 H h fc d )] + B. (A2) 

{fe*} 



7 i. 



on 



(A3) 



. . + 0(Xi,X 2 , 



, x d + Ax) 



+ 4>{x 1 ,x 2 , 



, x d - Ax)] 



By using Eq. I|A1(I in this last equation we obtain 

d 



cosh I A 



dxi 



b r =J2A(k*)J2cos(k % Ax)+B. (A4) 

i=l v {fc*} i=l 

As for n — 2, we note that 4 sinh 2 (y/2) sinh 2 (z/2) = [cosh(y) — l][cosh(z) — 1] and, in turn, cosh(y) cosh(;z) = 
i[cosh(y + z) + cosh(y — z)\. The latter combination leads to contributions that involve both forward and backward 
translations in different spatial directions. This is easily visualized by noting explicitly that 



cosh [ Ax— — 
' » dxi 



£cosh(Ax(A + A 



cosh Ax 



d 

dxi dxj 



d 



(A5) 



Notice that for the d cases where with i = j, the second term on the right hand side leaves the field at the original 
site r. The field at the original site is not represented by the ansatz assumption, and therefore we must subtract the 
d "spurious" terms produced by the ansatz state and add d times the field <f> r . This procedure leads to, 



cosh I A 



_d_ 

dxi 



d 
2' 



{fc*} 



cos(fcjAa 



-B. 



(A6) 



Note that we have taken advantage of the directional insensitivity of k* . 

Use of Eqs. I|A4(1 and l|A6() in Eq. 10 then leads to the following approximation for the term containing the Swift- 
Hohcnbcrg coupling operator: 



B{D\ — DHq), 



(A7) 



17 



where Di is given in Eq. (|13fl . 

Finally, since the summand in Eq. I|A7(I is independent of the direction of the k* , the sums simply give the number 
of terms n(fc*) in the sum, as given in Eq. IjlOjl (or the appropriate integral form), times the summand. Thus we 
finally arrive at the mean field approximation 

C4> r = Di [n(k*)A(k*) - 4> r ] + B{D X - Dk$). (A8) 

APPENDIX B: NORMALIZATION OF THE FOURIER TRANSFORM AND ORDER PARAMETER 

In this appendix we provide details on the relations between different relevant parameters used in the characteriza- 
tion of pattern formation. 

The Fourier transform of a field ip r and its inverse read respectively 



= Cj2^ lkT , (Bl) 

I" 

- C^ k e* r , (B2) 



k 

where C and C are normalization constants. Since 



^E e * k(r ~ r,) = ( B3 ) 

k 

^E^ (k - k0 - (B4) 

r 

the functional relation between C and C can be obtained by substituting Eq. I|B1|) into Eq. JB2}. One readily obtains 

d =cW- (B5) 

For simplicity we choose C = 1/N d , and therefore (7 = 1. 

Two parameters commonly used to characterize spatial patterns are the total power spectrum, S (fc), and the flux 
of convective heat, J, 

S(k) = ]T0 k 0- k , (B6) 
{k} 

J = jJsEtf. ( B7 ) 

r 

where the sum in (|B6|I runs over all modes of magnitude k. It is easy to check that the functional relation between 
these two quantities is simply 

J = Y,S(k), (B8) 

k 

where now the sum runs over the magnitude of the modes. 

For a pure spatial pattern of wave vector of magnitude k* where all spatial directions contribute in the same way 
we expect <fi T to be 

^ = E A ^ cos ( k ■ r ) • ( B9 ) 
{k*} 

Therefore, the Fourier transform of such a field is 

? k = ^(r)^(5 kjk .+5_ kjk ,). (BIO) 
{k*} 
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Consequently, the total power spectrum and the flux of this field are respectively, 



S(k) = -^>l 2 (r)^(5 k , k .+«5_ k , k .) 2 = n(fc*)^ 2 (fc*)^ ;fc .. (Bll) 
{k} {fe-} 

J = n(k*)A 2 (k*). (B12) 
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